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I 

SUMMARY  i 


This  report  summarizes  the  findings  of  an  experimental  study  supported  by  the  U.S. 
Army  Centrifuge  Research  Center  and  Engineer  Earthquake  Engineering  Research 
Program  (EQEN)  into  the  behavior  of  saturated  sands  under  high  initial  effective 
confining  stresses  subjected  to  strong  ground  shaking.  The  research  was  conducted 
using  the  Army  Centrifuge  at  the  U.S.  Army  Engineering  Research  and  Development 
Center  (ERDC),  located  in  Vicksburg,  MS.  A  large  dataset  of  the  response  of 
saturated  sand  to  dynamic  shaking  under  ‘level  groimd’  conditions  has  been  compiled 
and  a  series  of  verification  models  have  been  conducted.  Several  techniques  were 
used  to  investigate  the  response  of  deep  soil  sites  (in  excess  of  70m)  including 
surcharges,  lowered  water  table  and  higher  acceleration  (gravities).  The  most  effective 
approach  was  to  test  specimens  at  higher  gravity,  as  the  interpretation  of  data  from 
specimens  tested  using  a  heavy  surcharge  on  the  groimd  surface,  or  a  deep  water  table 
has  proved  difficult.  Many  deep  samples  showed  a  limit  to  the  excess  pore  pressure 
generated  during  shaking  at  a  level  much  less  than  100%  of  the  initial  vertical 
effective  stress.  A  range  of  explanations  has  been  proposed,  including  container 
effects,  saturation,  compression  of  the  soU,  redistribution  of  pore  pressure  and 
dynamic  response  of  the  surcharge.  The  verification  tests  eliminated  container  effects, 
and  parallel  research  by  others  has  shown  that  there  are  several  reasons  for  pore 
pressure  generation  at  depth  to  be  limited.  However,  as  postulated  in  earlier  stages  of 
this  research,  there  is  no  absolute  limit:  in  a  fuUy  saturated  specimen  of  broadly 
uniform  permeability,  liquefaction  was  observed  within  a  few  cycles  of  strong 
shaking  throughout  (to  an  equivalent  field  depth  of  around  65m).  A  second  key 
finding  indicates  that  dense  layers  overlying  loose  layers  may  still  be  readily  liquefied 
as  a  consequence  of  the  high  excess  pore  pressures  generated  below.  These  findings 
are  highly  significant  for  designers. 
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1.0  BACKGROUND 

This  research  contract  addressed  a  further  phase  in  the  completion  of  an 
experimental  and  analytical  research  programme  in  support  of  the  Earthquake 
Engineering  Research  Program  under  the  direction  of  the  USAE  ERDC  (formerly 
the  Waterways  Experiment  Station),  Vicksburg,  Mississippi.  The  phase  was 
entitled  Earthquake  Induced  Liquefaction  in  Deep  Sand  Deposits.  This  study  was  a 
continuation  of  earlier  research  under  Contract  Nos.  N68171-98-C-9014,  N68171- 
97-M-5710,  N68171-97-C-9012,  N68171-99-C-9021  and  N68171-00-M-5505. 

The  research  involved  further  interpretation  of  the  early  experimental  data,  design 
of  verification  experiments,  interpretation  and  analysis  of  new  experimental  data 
and  preparation  of  technical  papers  for  publication.  The  work  was  completed  in 
September  2001. 


2.0  INTRODUCTION 

The  current  state-of-practice  for  the  evaluation  of  liquefaction  potential  and  for 
remediation  design  and  analysis  depends  on  empirical  correlations  of  in-situ 
measurements  of  strength  versus  field  experience  of  liquefaction  at  shallow  depth 
and  laboratory  data  of  the  behavior  of  confined  elements  under  cyclic  loading. 
(Liquefaction  is  defined  here  to  mean  the  development  of  pore  pressure  equal  to 
100%  of  the  initial  vertical  effective  stress.)  This  approach  is  known  as  the 
“simplified  procedure”.  Opinions  vary  as  to  the  maximum  depth  in  the  field  at 
which  liquefaction  has  been  observed,  but  there  is  no  established  field  evidence 
fi-om  historic  earthquakes  of  liquefaction  at  depths  greater  than  a  few  tens  of 
meters.  The  NCEER  Workshop  in  1996  on  the  Evaluation  of  Liquefaction 
Resistance  of  Soils  noted  that  the  simplified  procedure  was  developed  fi’om 
evaluations  of  field  observations  and  field  and  laboratory  test  data,  Youd  and  Idriss 
(1997).  The  report  notes,  “These  data  were  collected  mostly  fi'om  sites  ...  at 
shallow  depths  (less  than  15m).  The  original  procedure  was  verified  for  and  is 
applicable  only  to  these  site  conditions”. 

Hence,  in  design  practice  the  assessment  of  liquefaction  imder  high  initial  effective 
confining  stress,  such  as  might  relate  to  the  foundations  of  large  earth  dams,  is 
based  on  the  extrapolation  of  observed  behavior  and  correlations  at  shallow 
depths.  In  practice,  the  behavior  of  saturated  soil  under  these  conditions  is  not  well 
understood.  Based  on  the  results  of  laboratory  tests,  researchers  have  postulated 
that  there  is  a  reduction  in  the  liquefaction  resistance  of  such  soil  compared  to 
shallow  depths.  This  reduction  is  accounted  for  in  standard  approaches  by  a  ratio 
known  as  K®,  a  “correction  factor”  developed  by  Seed  (1983)  in  the  simplified 
procedure.  This  strength  ratio  is  postulated  to  reduce  with  increasing  initial 
effective  confining  stress,  which  has  a  large  impact  potentially  on  the  extent  and 
method  (hence  cost)  of  remedial  construction  required  to  assure  adequate  seismic 
performance  of  large  dams.  For  example,  it  may  reduce  the  cyclic  shear  stress  ratio 
predicted  to  cause  liquefaction  in  a  soil  layer  under  a  typical  large  dam  (of  the 
order  of  30m  high)  to  about  50%  of  its  value  in  the  absence  of  the  dam.  The 
predicted  deformations  and  resistance  and  the  remedial  strength  required  are  a 
direct  function  of  the  reduced  shear  strength.  The  Ka  strength  reduction  therefore 
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has  a  strong  influence  on  the  decision  to  remediate  a  large  dam  and  on  the  costs  of 
that  remediation. 

The  fector  K<y  quantifies  the  curvature  in  the  cyclic  shear  strength  envelope  (cyclic 
shear  stress  required  to  cause  liquefaction  versus  confining  stress)  for  a  soil  as 
observed  in  lalwratory  tests  on  discrete  specimens.  Although  some  curvature  may 
be  expected,  such  large  reductions  in  cyclic  shear  strength  ratios  are  counter¬ 
intuitive.  It  is  generally  accepted  that  increased  confining  stress  should  broadly 
improve  the  capacity  of  a  soil  to  resist  applied  loads,  not  reduce  it.  Clearly  the 
volume  of  soil  that  requires  to  be  treated  and  the  difficulty  and  expense  of  that 
treatment  are  highly  dependent  on  an  accurate  assessment  of  the  potential  for  and 
consequences  of  liquefaction. 

There  is  therefore  strong  motivation  for  owners  of  large  dams  to  investigate  the 
behavior  of  saturated  sands  subject  to  strong  ground  shaking  under  high  initial 
effective  confining  stresses.  In  the  absence  of  field  data,  the  use  of  a  centrifuge 
was  considered  to  be  the  only  practical  option  to  realistically  represent  a  deep  soil 
deposit  subjected  to  earthqu^e  shaking.  The  experiments  reported  herein  relate 
exclusively  to  level  ground  initial  stress  conditions  in  a  single  layer  (loose)  or  two- 
layer  (dense  over  loose)  deposit  of  clean,  fine  sand. 

Liquefaction  of  a  level  sand  bed  has  previously  been  the  subject  of  research  by  many 
others.  Experiments  using  Nevada  sand  were  conducted  at  many  centrifuge  centers 
under  the  VELACS  project,  Aralanandan  and  Scott  (1994).  The  objective  of  the 
ERDC  study  is  to  investigate  the  onset  of  liquefaction  imder  much  higher  initial 
effective  overburden  stresses  than  used  in  the  VELACS  experiments.  Nevada  sand, 
which  was  used  in  the  VELACS  research  and  has  been  extensively  studied  in  the 
laboratory,  was  adopted  for  the  series  of  e^qjeriments  under  the  ERDC  centrifuge  test 
program. 


3.0  THE  ERDC  CENTRIFUGE  FACILITY 

The  design  specification  for  the  ERDC  centrifuge  followed  a  review  of  the  Army’s 
research  needs  and  a  study  of  the  available  academic  facilities,  Ledbetter  (1991). 
Many  of  the  field  problems  with  which  the  U.S.  Army  Corps  of  Engineers  is 
concerned  are  physically  large,  such  as  earth  dams,  locks  and  river  control 
structures,  environmental  problems  and  military  research.  It  was  determined  that  a 
new  facility  with  a  high  payload  capacity  and  high  g  capability  was  required  to 
meet  future  Army  needs.  A  large  beam  centrifuge  was  commissioned.  Figure  1, 
termed  the  Acutronic  684,  based  on  the  French  designed  Acutronic  661,  665  and 
680  series  of  geotechnical  centrifuges,  Ledbetter  et  al.  (1994a).  The  capacity  of  the 
ERDC  centrifuge  is  a  payload  of  8  tonnes  at  up  to  143  g,  reducing  to  2  tonnes  at 
350  g,  with  a  platform  area  of  1.3m  square.  This  high  capacity  enables  field 
problems  of  the  order  of  up  to  300m  in  breadth,  300m  in  depth,  and  1 000m  in 
length  to  be  simulated  under  a  wide  variety  of  loading  conditions.  The  facility  is 
equipped  with  a  large  range  of  equipment  and  appurtenances. 
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Figure  1.  The  WES  centrifuge 

The  research  approach  for  this  high  confining  stress  liquefaction  study  uses  the 
large  capacity  of  the  WES  centrifuge  to  investigate  the  generation  of  excess  pore 
pressure  and  liquefaction  under  conditions  that  much  more  closely  resemble  those 
at  depth  in  the  field,  Ledbetter  et  al.  (1999).  To  do  this  required  the  design  and 
construction  of  a  large  dynamic  actuator.  The  history,  design  and  characteristics 
of  the  actuator  and  specimen  container  have  been  discussed  in  a  previous 
Technical  Report,  Steedman  (2000)  and  are  not  repeated  here. 

4.0  RESEARCH  PROGRAM 

Table  1  summarises  the  experiments  conducted  during  1998  and  2001.  The  models 
are  grouped  in  series,  where  each  series  corresponds  to  a  different  target  range  of 
vertical  effective  overburden  stress  in  the  loose  layer. 


Model 

series 

Models  in  series 

Effective 
overburden  stress 
near  bottom 

Equivalent  field 
depth  (approx) 

Depth  of 
specimen 

Notes  (all  specimens  constructed  from 
Nevada  Sand,  tested  at  50g  unless  specified) 

2 

a,  b,  c,  d,  e,  f 

llOKPad  tsf) 

15  m 

300  mm 

Fully  submerged 

3 

a,  b,  c,  d,  e,  f 

210KPa(2tsf) 

26  m 

525  mm 

Fully  submerged  (3f  at  lOOg) 

4 

a,  b,  c,  d 

320-540KPa 
(3  -  5  tsf) 

26  -  40  m 

525  mm 

Lowered  phreatic  surface  or  surcharge 
(weights) 

e,f,g,h,  i,j 

420- 1040  KPa 
(3.9 -9.7  tsf) 

30-60m 

525  mm 

Surcharge  (weights)  and  fully  submerged, 
tested  up  to  125g 

k 

200 -510  KPa 
(1.9 -4.8  tsf) 

25  -  39  m 

525  mm 

No  surcharge,  fiilly  submerged,  tested  up  to 
125g 

5 

a,  b,  c,  d,  e 

750- 1070  Kpa 
(7 -10  tsf) 

54  -  63  m 

525  mm 

Surcharge  (weights),  fully  submerged 

5  g 

740  KPa  (6.9  tsf) 

41  m 

525  mm 

Surcharge  (weights),  fully  submerged  (5g 
reused  Model  3e) 

h 

610ICPa(5.7  tsf) 

41  m 

525  mm 

No  surcharge,  fully  submerged,  tested  up  to 

125g 

i 

730  KPa  (6.8  tsO 

41  m 

460  mm 

No  surcharge,  lowered  phreatic  sur&ce, 
tested  at  125g 

Table  1.  Outline  summary  of  centrifuge  model  test  program 
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Overall 
depth  (mm) 

Relative  Density  (Dr)* 

Ov’  near  base  of 
specimen  (tsf) 

OCR 

Number  of 
earthquakes 

Comments  (all  models  constructed  from  Nevada 
Sand  and  tested  at  50g  unless  specified) 

2a 

300 

44%  loose,  83%  dense 

1 

1 

3 

Fully  submerged 

2b 

300 

50%  loose,  75%  dense 

1 

1 

2 

Fully  submerged 

2c 

300 

49%  loose,  74%  dense 

1 

1 

5 

Fully  submerged 

2d 

300 

50%  loose,  75%  dense 

1 

1 

4 

Fully  submerged 

2e 

300 

49%  loose,  73%  dense 

1 

2.5 

4 

Fully  submerged 

2f 

300 

50%  loose,  75%  dense 

1 

2.5 

4 

Fully  submerged 

3a 

525 

34%  loose,  73%  dense 

2 

1 

2 

Fully  submerged 

3b 

525 

49%  loose,  77%  dense 

2 

1 

3 

Fully  submerged 

3c 

525 

49%  loose,  79%  dense 

2 

1 

3 

Fully  submerged 

3d 

525 

54%  loose,  80%  dense 

2 

2.5 

4 

Fully  submerged 

3e 

525 

50%  loose,  75%  dense 

2 

1 

1 

Fully  submerged 

3f 

262 

50%  loose,  75%  dense 

2 

1 

1 

182  mm  dense  layer  overlying  80  mm  loose 
layer,  fully  submerged,  no  surcharge,  lOOg 

4a 

525 

49%  loose,  80%  dense 

3 

1 

4 

Saturated  to  top  of  loose  layer  only. 

4b 

525 

56%  loose,  74%  dense 

3 

2.5 

4 

Saturated  to  top  of  loose  layer  only. 

4c 

525 

50%  loose,  75%  dense 

4.7 

1 

4 

Fully  submerged,  surcharge  (weights) 

4d 

525 

50%  loose,  68%  dense 

4.7 

2.5 

4 

Fully  submerged,  surcharge  (weights) 

4e 

525 

47%  uniform 

3.9 -7.8 

1 

2,1,1 

Fully  submerged,  surcharge  (weights).  Shaking  at 
50,  80,  lOOg 

4f 

525 

55%  uniform 

3.9 -9.7 

1 

1,1, 2,1 

Fully  submerged,  surcharge  (weights).  Shaking  at 
50,  80, 100, 125g 

4g 

525 

50%  uniform 

3.9  -  9.7 

1 

1,1,  1,2 

Fully  submerged,  surcharge  (weights).  Shaking  at 
50,  80, 100,  I25g 

4h 

525 

50%  uniform 

3.9  -  9.7 

1 

1,1, 1,1 

Fully  submerged,  surcharge  (weights).  Shaking  at 
50,  80, 100,  125g 

4i 

525 

50%  uniform 

3.9  -  9.7 

1 

1,1,  1,1 

Fully  submerged,  surcharge  (weights).  Shaking  at 
50,  80, 100,  125g 

4j 

525 

50%  uniform 

3.9 -9.7 

1 

1,1,  1,2 

Fully  submerged,  surcharge  (weights).  Shaking  at 
50,  80, 100, 125g. 

4k 

525 

50%  uniform 

1.9-4.8 

1 

1,1, 1,2 

Fully  submerged,  no  surcharge.  Shaking  at  50, 

80, 100,  125g. 

41 

525 

Densities  uncertain 

4.8 

1 

1 

Fully  submerged,  surcharge  (weights) 

5a 

525 

5 1%  loose,  72%  dense 

7.4 

1 

4 

Fully  submerged,  surcharge  (weights) 

5b 

525 

49%  loose,  76%  dense 

7.4 

2.5 

4 

Fully  submerged,  surcharge  (weights) 

5c 

525 

52%  loose,  75%  dense 

9.2 

1 

3 

Fully  submerged,  surcharge  (weights) 

5d 

525 

57%  loose,  80%  dense 

9.2 

1 

1 

Fully  submerged,  surcharge  (weights) 

5e 

525 

50%  uniform 

8.4 

1 

7 

Fully  submerged,  surcharge  (weights) 

5f 

525 

50%  loose,  75%  dense 

6.9 

1 

1 

Fully  submerged,  surcharge  (weights) 

5g 

525 

50%  loose,  75%  dense 

6.9 

1 

1 

Fully  submerged,  surcharge  (same  model  as  3e 
with  surcharge  added) 

5h 

525 

49%  loose,  74%  dense 

5.7 

1 

1 

Fully  submerged,  no  surcharge,  125g 

5i 

460 

50%  loose,  75%  dense 

6.8 

1 

1 

396  mm  dense  layer,  saturated  to  depth  of 

250mm  only,  overlying  64  mm  loose  layer,  no 
surcharge,  125g 

♦Relative  density  of  the  portion  of  the  model  termed  Moose’,  generally  refers  to  a  160  mm  thick  layer  placed  at  the  bottom 
of  the  model.  For  the  ‘model  of  the  model’  experiments,  this  loose  layer  was  reduced  in  depth.  Relative  density  of  the 
portion  of  the  model  termed  ‘dense’,  refers  to  that  portion  of  the  model  above  the  ‘loose’  layer  (ranged  in  thickness  from 
140  to  365  mm).  Model  41  was  understood  to  have  b^n  built  to  the  standard  specification,  i.e.  75%  dense  over  50%  loose. 

Table  2.  Summary  description  of  centrifuge  model  test  program 
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Generally,  specimens  were  up  to  525mm  deep  with  the  bottom  160mm  placed  at 
around  50%  Relative  Density  (RD)  and  the  upper  portion  at  around  75%  RD.  The 
‘early’  models  were  shaken  at  50g;  later  models  were  shaken  at  up  to  125g.  During 
the  verification  phase  modelling  of  models  was  carried  out  at  lOOg  and  125g. 

Some  models  were  overconsolidated  by  a  factor  of  2.5  prior  to  shaking  (achieved 
by  running  the  centrifuge  up  to  125g  and  then  dropping  back  to  50g).  The  size  of 
the  dataset  and  number  of  repetitions  provided  the  opportunity  to  characterise  the 
data  on  the  basis  of  the  quality  of  the  output.  This  is  seen  in  Appendix  A,  where 
the  summary  data  for  each  transducer  and  each  earthquake  analysed  is  presented. 

A  quality  rating  of  low/medium/high  and  *  was  assigned  to  each  time  history 
based  on  the  consistency  and  reliability  of  the  data  and  the  quahty  of  the 
experiment  itself  For  example,  higher  ‘quality  ratings’  would  be  assigned  to 
transducers  which  showed  internal  consistency  in  terms  of  static  pore  pressures, 
physical  location  in  the  model,  response  during  shaking  etc. 

The  Nevada  sand  used  in  the  models  was  characterized  by  standard  laboratory 
tests  to  determine  parameters  such  as  dry  density  and  gradation.  Table  3  presents 
key  material  parameters  for  this  sand. 


Specific  gravity 
Maximum  void  ratio 
Minimum  void  ratio 

Dso 

Dio 


2.64 

0.757  (density  93.8  pcQ 
0.516  (density  108.7  pcf) 
0.18  mm  (approx) 

0.11  mm  (approx) 


Table  3.  Nevada  Sand  (parameters  as  measured) 

The  pore  fluid  comprised  a  mixture  of  glycerine  and  water,  80%  by  weight  for 
experiments  conducted  at  50g.  Measurements  of  the  viscosity  of  glycerine-water 
mixes  at  a  range  of  temperatures  and  proportions  show  that  the  viscosity  is 
sensitive  to  both  parameters. 

The  density  of  a  glycerine-water  mix  was  calculated  from: 


pm  =  pg(mg  -I-  mw)/(  mg  +  pg  mw) 


where  pm  is  the  density  of  the  mix,  pg  is  the  density  of  glycerine,  nig  is  the  mass  of 
glycerine,  and  mw  is  the  mass  of  water.  Table  4.  summarizes  the  properties  of  the 
glycerine-water  solution  used  as  the  pore  fluid. 

Density  1200  kg/m^ 

Viscosity  50  cs 

Specific  Gravity  1.26 

Composition  80%  glycerine-water  mix  (by  weight) 

Table  4.  Parameters  for  pore  fluid  (as  measured) 
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The  models  were  poured  dry  from  a  hopper  and  saturated  under  vacuum,  or  slowly 
under  gravity.  Instrumentation  was  placed  in  the  model  as  it  was  being 
constructed. 


The  typical  arrangement  of  the  model  specimens  for  the  different  test  series  is 
shown  in  Figure  2.  Instrumentation  was  positioned  through  the  depth  of  the 
models,  and  comprised  pore  pressure  transducers  and  accelerometers.  The 
location  of  the  instrumentation  for  Model  5b  is  shown  in  Figure  3. 


Typical  cross  section  through  Typical  cross  section  through  Model  series  3 


Typical  cross  section  through  Models  4a,  4b  Typical  cross  section  through  Models  4c,  4d  and  5 


Figure  2.  Cross  sections  through  the  different  model  configurations 
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Figure  3.  Instrumentation  for  Model  5b 


5.0  RESULTS 

Earlier  reports,  Steedman  (1999),  (2000)  have  presented  analysis  of  the  data  using 
a  range  of  approaches.  In  this  report,  the  data  from  the  verification  models  is 
summarized  alongside  the  earlier  experimental  results,  and  in  the  same  format. 

Appendix  A  presents  the  dataset  of  limiting  excess  pore  pressures,  by  model  and 
by  transducer.  Comments  are  added  concerning  the  quality  of  each  record,  and  the 
experiment  in  general. 

From  this  set  a  number  of  time  histories  were  selected  as  of  particularly  high 
quality  for  possible  future  interpretation.  This  selection  was  made  in  part  based  on 
the  internal  consistency  of  the  experimental  results,  and  in  part  using  subjective 
judgment  based  on  past  experience.  These  data  are  presented  in  Figures  4a  and  4b, 
and  5a  and  5b,  below.  In  each  case,  the  model  reference  number  and  earthquake  is 
given,  which  can  be  checked  against  the  commentary  in  Appendix  A.  Figures  4a 
and  5a  present  the  pore  pressure  records;  Figures  4b  and  5b  present  the  (matching) 
base  input  acceleration  time  history  from  the  same  model  test. 

A  summary  of  selected  data  and  recommendations  for  future  centrifuge  testing  of 
level  ground  conditions  under  high  effective  confining  stress  has  been  prepared, 
Sharp  and  Steedman  (2002). 
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Figure  4a.  Selected  pore  pressure  time  histories  (I) 
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Figure  5a.  Selected  pore  pressure  time  histories  (II) 
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6.0  DETERMINATION  OF  SHEAR  STRESS-STRAIN  HISTORIES 

All  of  the  experiments  undertaken  in  this  study  addressed  level  ground  conditions, 
with  no  initial  static  shear  stress.  The  only  shear  stress  and  subsequent  shear  strain 
was  that  imposed  by  the  loading.  Level  ground  liquefaction  is  therefore  a  special 
case  of  the  cyclic  mobility  phenomenon  that  occurs  when  the  static  shear  stress  is 
less  than  the  shear  strength  of  the  liquefied  soil.  Level  groxmd  liquefaction  failures 
are  caused  by  the  upward  flow  of  water  that  occurs  when  seismically  induced 
excess  pore  pressures  dissipate.  Depending  on  the  length  of  time  required  to  reach 
hydraulic  equilibrium,  level  ground  hquefaction  failure  may  occur  well  after 
ground  shaking  has  ceased  and  this  was  seen  in  a  number  of  the  experiments. 
References  in  Appendix  A  to  transducers  ‘still  climbing  at  the  end  of  shaking’  are 
evidence  of  pore  pressure  redistribution  affecting  the  pore  pressure  generation 
caused  by  shaking. 

A  critical  aspect  in  the  assessment  of  these  model  experiments  was  to  compute  the 
stress  strain  time  histories  for  comparison  with  laboratory  element  and  field 
experience.  This  was  done  in  collaboration  with  researchers  at  WES  and  RPI. 

The  method  adopted  to  compute  the  stress  and  strain  time  histories  was  known  as 
System  Identification  or  SI.  The  method  is  described  in  a  paper  by  Zeghal  et.  al. 
(1995),  where  the  authors  proposed  that  accelerometer  and  pore  pressure  time 
histories  could  offer  a  direct  and  effective  means  of  evaluating  seismic  soil 
properties,  in  particular,  a  second  order  accurate  estimation  of  shear  stress  and 
shear  strain.  This  SI  technique  was  utilized  for  analysis  of  the  earthquake  data. 

Stress  strain  figures  are  shown  here  for  Models  4k  (Figure  6)  and  5b  (Figure  7), 
windowed  during  similar  stages  of  shaking,  Steedman  and  Sharp  (2001). 

The  left  column  of  plots  are  the  traditional  stress-strain  loops  computed  for  the 
early  time  data,  as  the  excess  pore  pressures  are  developing  strongly.  This  may  be 
seen  by  comparison  with  the  time  history  of  vertical  effective  stress  below  the 
stress  strain  figures.  For  Model  4k  this  shows  the  largest  amount  of  straining  near 
the  surface  at  nearly  1%,  decreasing  to  a  strain  of  0.6%  at  a  depth  of  24  m.  The 
remaining  columns  of  figures  show  how  the  stress  strain  loops  are  altered  by  the 
development  of  excess  pore  pressure.  In  the  case  of  Model  4k,  the  specimen 
reaches  close  to  100%  excess  pore  pressure,  and  the  upper  layers  become  isolated 
jfrom  the  straming  at  the  base.  In  the  case  of  Model  5b,  the  excess  pore  pressure  is 
limited  and  the  response  of  the  soil  column  becomes  nearly  linear.  The  marked 
contrast  in  response  between  the  two  specimens  confirm  the  observation  that  the 
stabilization  of  the  excess  pore  pressure  development  in  the  deeper  specimens  is  a 
genuine  phenomenon  associated  with  a  stable  cycling  of  applied  load  at  a  level 
significantly  below  100%  excess  pore  pressure,  and  with  the  soil  retaining 
considerable  reserves  of  strength  and  stif&iess. 

One  explanation  for  the  limiting  level  of  excess  pore  pressure  may  be  associated 
with  upward  drainage  of  the  excess  pore  pressures. 

Appendix  B  presents  isochrones  of  stress  and  strain  levels  from  each  experiment, 
based  on  this  method. 
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A  very  large  dataset  of  the  behavior  of  Nevada  sand  subject  to  shaking  under  level 
ground  conditions  has  been  accumulated  as  a  result  of  this  work.  Review  of  the 
stress-strain  data  alongside  the  time  histories  of  excess  pore  pressure  provided 
considerable  additional  insight  into  the  explanation  of  the  pore  pressure  response 
in  many  cases.  The  verification  models  that  were  tested  during  2000-2001  were 
successful  in  validating  the  performance  of  the  specimen  chamber  and  shaker  in 
many  respects.  The  verification  models  showed  that  the  use  of  surcharges  may 
indeed  be  having  an  effect  on  the  response  of  the  soil  column,  although  exactly 
how  the  surcharge  affects  the  soil  column  has  not  yet  been  conclusively 
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determined;  many  other  factors,  including  soil  density,  compression  and  saturation 
may  also  be  influencing  the  ‘true’  response  under  cyclic  loading  at  high  effective 
confining  stress.  It  is  clear  that  the  most  effective  method  of  investigating  the 
influence  of  very  high  effective  confining  stress  is  to  use  the  centrifuge  at  high  g, 
with  or  without  dry  sand  surcharging,  as  for  example  was  achieved  with  the  final 
experiments  in  the  series.  Models  5h  and  5i. 
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Figure  8.  Liquefaction  at  62m  equivalent  depth.  Model  5h 

Figure  8  presents  data  from  PPT35,  Model  5h,  which  shows  100%  excess  pore 
pressure  reached  after  a  few  cycles  of  strong  shaking  in  a  fully  saturated  deep 
specimen  without  surcharge,  at  125g.  The  shear  stress  ratio  time  history  shown 
here  was  calculated  using  the  simplified  method,  Steedman  (2000)  and  may  be 
compared  to  the  isochrones  of  stress  ratio  calculated  using  the  SI  method, 
Appendix  B.  Appendix  B  also  shows  the  long  term  response  after  the  end  of 
shaking,  where  consolidation  effects  and  pore  pressure  redistribution  last  for  a 
very  long  period  after  the  end  of  the  earthquake. 
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7.0  CONCLUSIONS 

1 .  A  large  data-set  of  the  behavior  of  loose  saturated  sands  under  high  initial 
effective  confining  stresses  and  subject  to  earthquake-like  shaking  has  been 
collected  during  an  extensive  experimental  program  on  the  ERDC 
Centrifuge,  Vicksburg  MS. 

2.  Early  findings  from  the  centrifuge  experiments  showed  that  under  high 
effective  confining  stress,  the  potential  to  generate  excess  pore  pressures 
was  limited,  and  this  may  be  of  considerable  importance  if  verified. 

3 .  Detailed  investigation  into  the  early  findings  has  provided  a  range  of 
explanations  for  the  limiting  values,  including  experimental  problems  such 
as  specimen  saturation,  chamber  effects  or  the  distorting  effects  of  large 
surcharges  on  the  soil  surface.  Of  these  possible  effects,  detailed 
investigation  of  the  response  of  the  chamber  effectively  eliminated 
movement,  or  distortion  of  the  ESB  container  as  an  explanation. 

4.  The  effect  of  surcharging  was  investigated  by  comparison  with  changing  g 
level  and  using  different  methods  of  surcharging,  or  no  surcharges  at  all. 
The  use  of  dry  sand  appeared  preferable  to  lead  or  steel  weights;  this 
method  was  used  in  early  experiments  as  well  as  in  later  models. 

5.  In  the  final  experiments,  which  were  carried  out  at  high  g  without  the  use 
of  surcharges,  liquefaction  occurred  throughout  the  entire  depth  within  a 
few  cycles.  This  experiment  confirmed  that  with  sufficient  input  shaking, 
large  excess  pore  pressures  are  generated  in  a  saturated  soil  column  of 
broadly  uniform  permeability. 

6.  A  key  finding  of  the  study  has  been  to  recognise  the  importance  of  pore 
pressure  redistribution  in  deep  soil  columns,  which  can  cause  an  ‘up- 
welling’  of  excess  pore  pressure  from  lower  layers  leading  to  liquefaction 
of  upper  layers  some  time  after  shaking  has  ceased.  This  is  ein  effect  not 
widely  considered  in  standard  design  methods. 
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APPENDIX  A 

EARTHQUAKE  DATASET 

Data  of  limiting  pore  pressures.  Models  and  cycles  (1998-2001) 
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APPENDIX  B 

SUMMARY  OF  DATA 

Time  histories  of  acceleration,  excess  pore  pressure  and  model  configurations 

Stress  and  strain  histories 
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Model  3e  Stress  and  strain  isochrones 
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Model  4k  Stress  and  strain  isochrones 
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